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We establish a technique to find the states with most robust entanglement in dissipative quantum 
systems and explicitly construct those state for various environments. 
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The potential to form entangled states is one of the 
central distinctions between quantum objects and their 
classical counterparts. Therefore, the stability properties 
of entanglement teaches us a lot about the emergence of 
classicality in quantum systems of growing size. Besides 
that - or also, even more importantly - entangled states 
are the central building block for many promising appli- 
cations [l[ in quantum information theory, so that there 
is a considerable interest not only in the preparation of 
highly entangled states, but also in the preservation of 
entanglement over sufficiently long times that allow to 
execute a quantum algorithm or different tasks. 

Generally, entanglement decays due to environment 
coupling [2, 0, HI- However, the decay of entanglement 
does not necessarily follow that of the density matrices 
coherences. In particular, since entanglement is not a lin- 
ear function of the underlying quantum state, there can 
be states whose entanglement is significantly more robust 
then that of other states. 

Here, we seek those states whose entanglement is most 
robust in a given situation of environment coupling. That 
is, given an entanglement measure E, we look for those 
states, for which the temporal increment E is maximal. 

Such tasks require an entanglement measures that al- 
lows a simple evaluation or estimation. Entanglement 
measures, however, are rather intricate to evaluate since 
most of them rely on a mathematical optimization pro- 
cedure that can be solved without extensive numerical 
optimizations only in exceptional cases. What we con- 
sider here, is not a measure E(g) itself, but rather its 
time derivative E(g,g), and the technical difficulties to 
evaluate the letter quantity are certainly not smaller than 
for the former. 

Therefore, we make use of a purely algebraic estimate 
of concurrence [f| in terms of a bilinear functional of the 
density matrix g. The concurrence of a pure state "J) can 
be defined via a spin-flip operation for two-level systems 
5], or in terms of a linearized entropy of the reduced 
density matrix of one of the subsystems [H, or, also, via 
the expectation value of a suitably chosen operator with 
respect to the duplicate quantum state \^} ® |\?) Q for 
systems of arbitrary finite dimension. The extension to 
mixed states can be performed with the help of a convex 
roof 0] c(g) = inf ^2iPi\^i)(^i\, where the infimum is 
to be found among all pure state decompositions g — 
YliPi\^i)(^i\ of the state g. Whereas this infimum can 
be found algebraically only for small systems 0] or some 



special states [10| . it can be bounded from below by 

c 2 (g) =r(g) >Trg®gV , (1) 

for general states [ll|, where V = P- <£> P_ - 1/2(P- ® 
P + + P + (S>P-) is defined in terms of the projectors P T on 
the antisymmetric and symmetric components of the du- 
plicated Hilbert spaces of either subsystem, and t = c 2 is 
called tangle. This is exact for pure states, and provides 
a very good approximation for weakly mixed states, what 
makes it an ideal tool for the following investigations in 
which we will consider the decay of initially pure states 
into mixed ones. With this bound, we can approximate 
the temporal increment of tangle as 



T(Q,Q)~2Tr(g® V) , 



(2) 



what allows to maximize f over pure initial states. Typ- 
ically, weakly entangled states show more robust entan- 
glement than highly entangled ones, and initially sepa- 
rable states have a vanishing temporal increment under 
coupling to local environments. In order to ensure that 
our optimizations will not simply yield separable states 
as those with most robust entanglement, we will fix the 
initial tangle t(\I/) = tq, and seek those states with given 
tq that maximize the temporal increment of tangle 
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max(f (*)|t(*) = r ) 



(3) 



Environment coupling will be described in the follow- 
ing in terms of a Master equation [12L 1 1 31 ] g = 2; C ai g, 
with a Lindbladian C a that acts like 



C a g = (2<jgo-' — a'ag — ga'a 



(4) 



For coupling operators a we consider the case of local en- 
vironments, where any such quantities acts nontrivially 
only on a single subsystem, since this corresponds to the 
typical situation in which the subsystems are macroscop- 
ically separated. The two iconic models for evolutions of 
quantum states in quantum information processing are 
the the 'amplitude damping' channel that corresponds 
to a = o~- = (o x — ia y )/2, and the 'phase damping' 
channel that corresponds to a — o~ z [14| in terms of the 
Pauli spin operators. Here, we consider generalizations 
of these situations to quantum systems with more than 
just two levels. Typical generalizations to higher dimen- 
sional systems are often derived for the harmonic oscil- 
lator, where the generalization of <r_ is the annihilation 
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FIG. 1: Schmidt coefficients of the states with the most robust entanglement against spontaneous decay as function of tangle 
r. The different insets correspond to different values of the parameter q that characterizes the decay mechanism in Eq. ([SJ. 
The Schmidt coefficients are distinguished by their dashing with increasing length of dashing from Ao (shortest dashing) to A3 
(solid line). 



operator a, and the generalization of a z is the number 
operator n = a^a. This, however, is a rather specific 
situation that can not describe general experimental sit- 
uations sufficiently well. For example, there could be 
several excited states that have comparable life time, or 
rather different ones that are incompatible with those re- 
sulting from the harmonic oscillator operators. We there- 
fore take the generalization 

^dc(9) = E( i + 1 ) 9 l^l ' ( 5 ) 

i 

of the harmonic oscillator case to describe decoherence. 
Here, there is the variable parameter '5' that allows to 
change the properties of the dephasing mechanism. The 
harmonic oscillator case is recovered for 5 = 1, but chang- 
ing the value of q allows to adjust the decay rates for the 
coherences between different levels. As generalization of 
spontaneous decay, we take 

a ad ( 9 ) = E(z + l) 9 |0)(z + l| , (6) 

i 

that is the case where all excited levels decay to the 
ground state with adjustable decay rates. This is differ- 
ent to the harmonic oscillator case, where excited states 
decay to the ground state via all lower lying states, but 
rather resembles the situation of an atom or ion with 
several excited states the decay directly to the ground 
state. 

Obviously, also these models can not provide an ex- 
haustive description of very general dissipative dynamics. 
However, these two situations will provide a good insight 
in the interplay of stability of entanglement and environ- 
ment coupling, and situations with specific differences as 
compared to Eqs. ([5]) and © can easily be investigated 
along the lines that we present in the following. 

Like any entanglement measure, the tangle is invariant 
under local unitaries, and given this invariance, any state 
can be represented in its Schmidt decomposition [151 ] 

$^VA^i>®l&>, (7) 

i 



where all entanglement properties are described in terms 
of the Schmidt coefficients Xi, and and {|0»}} are 

local bases for the individual subsystems. However, the 
environment coupling typically breaks this invariance for 
f, so that the present maximization can not be restricted 
to states that are of Schmidt form (Eq. ([7])). Never- 
theless, most situations of environment coupling define 
a special basis, the so-called pointer basis that consists 
of states that are most robust under environment cou- 
pling. Performing numerical optimizations over general 
states we found that typically states with such pointer 
bases as Schmidt bases have the most robust entangle- 
ment. We will therefore take exactly those states, op- 
timize over their Schmidt coefficients, and convince our- 
selves at the end with comparison to numerical optimiza- 
tions over general states that this restriction still allows 
to find the states with the most robust entanglement. 

Given an initial state in Schmidt decomposition, the 
most robust entanglement is determined by the solution 
of polynomial equations 

E ai ^xP = < andr(*)=r (8) 

i 

where the prefactors oti need to satisfy X^^TTa^ = 
to ensure that the variation over the Xi is performed over 
states with a constant tangle. Such sets of polynomial 
equations can very reliably and efficiently be solved with 
the help of Groebner bases [l(| . The risk of finding local 
instead of global minima that always affects numerical 
optimization is basically dispelled. 

Fig.Q]shows the Schmidt coefficients of the states with 
the most robust tangle against environment coupling as 
given in Eq. $Q together with Eq. (01 as function of tan- 
gle r for different values of the parameter q. For tangles 
larger than 1 , necessarily three Schmidt coefficients need 
to be positive, and values larger than 4/3 require at least 
four positive Schmidt coefficients. The right-most inset 
of Fig. [T] shows the a situation with significantly different 
decay times for the individual Schmidt-basis states. A 
third and fourth basis state are occupied only for a value 
of r above the threshold values r = 1, and r = 4/3. This 
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FIG. 2: Schmidt coefficients of the states with most robust entanglement against decoherence analogous to Fig.[T] The different 
insets correspond to different values of q defined in Eq. Jjj). 



is different in the case of comparable life times of the in- 
dividual basis states as depicted in the left- most inset. 
Here, a third and fourth level are occupied already for 
significantly smaller values of r. That is, even if a given 
value of tangle can be realized also with only two occu- 
pied levels, it is actually favorable to occupy more levels, 
and thereby enhance the stability of entanglement. This 
enhancement can be attributed to a larger occupation of 
the stable ground state, as it can be seen in Fig. [T] when- 
ever an additional level gets occupied, there is a kink in 
the largest Schmidt coefficient that shows how the de- 
crease of the population of the stable levels is slowed 
down. 

The situation is quite different in the case of decoher- 
ence that is displayed in Fig. O For values of q that are 
smaller than one, it is indeed favorable to occupy the 
levels 1 2} and |3), since the decoherence times for super- 
positions of the states |2) g> |2) and |3) <g> |3) are longest. 
If q is larger than one, the situation changes and super- 
positions of the states |0) ® |0) and |1) ® |1) are most 
stable so that there is a qualitative change between the 
two insets that correspond to the values q — 24/25 and 
q = 26/25. A third and fourth level is occupied only 
above the threshold values r = 1 and r — 4/3 for all 
values of q. However, the distribution of Schmidt coeffi- 
cients at a given tangle r > 1, rather strongly depends on 
q: whereas the large Schmidt coefficients are of compara- 
ble size while there is single small Schmidt coefficient for 
q -C 1 and for q 3> 1, it is more favorable to have more 
broadly distributed Schmidt coefficients for q ~ 1. 

So far, we have been assuming that the optimal 
states have their Schmidt bases given in terms of the 
environment-induced pointer bases. Fig. [3] shows a com- 
parison of the algebraic solutions obtained under this as- 
sumption with numerical solutions, where the optimiza- 
tion has been performed over general initial states. The 
left inset displays the logarithm log(— r) of the temporal 
increment f of tangle for different values of q. Appar- 
ently, all numerically obtained data points lie above the 
algebraic solutions that are shown as lines. This implies 
that no improvement in the optimization can be obtained 
dropping the assumption on the Schmidt bases. Besides 
that, there are also several data points that lie signifi- 
cantly above the algebraic solutions. These points cor- 
respond to local maxima that pose a serious challenge 
to numerical optimizations, that, however, are hardly a 
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FIG. 3: Comparison of numerically optimized temporal in- 
crement of tangle and optimized f with pointer bases as 
Schmidt bases: the left inset corresponds to the situation of 
dephasing and the right inset to spontaneous decay. The cir- 
cles, squares and triangles display the numerical solutions for 
q = 1/4,1/2,3/4,1,5/4 and 3/2 (from bottom to top), and 
the lines show the algebraically obtained solutions as function 
of the parameter q defined in Eqs. (|6]) and |[5j. 



problem for our algebraic solutions. The right inset of 
Fig. [3] shows the analogous situation for decay. Here, 
the numerical optimization is significantly more reliable 
than in the case of decoherence, and there is only one 
data point that corresponds to a local maximum. Thus, 
the conjecture that Schmidt and pointer bases coincide 
is justified by the data shown. 

Finally, let us compare the optimal solutions that have 
been obtained with the approximation Eq. (J5]) with the 
time-dependence of state that is not subject to this ap- 
proximation. Here, we use the bounds [l?], HH that are 
known to provide a very accurate estimate of concur- 
rence, in particular for weakly mixed states, this is ex- 
actly the situation that we are dealing with here. The 
black line in the left plot of Fig. 2] shows the time evo- 
lution of r for the state with initial tangle tq = 7/3 for 
decoherence with q = 1. The grey lines show the tan- 
gle for 1000 randomly chosen states with the same initial 
tangle. The tangle of the optimized case clearly decays 
significantly slower than that of the other states, and 
the optimized states carry sizable entanglement at times 
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FIG. 4: Time dependence of tangle for the obtained most ro- 
bust states (black line) an 1000 randomly chosen initial states 
with a given initial tangle To- The left plot corresponds to de- 
coherence (Eq. ([5])) with to = 7/3, the right plot to decay 
(Eq. @) with ro = 8/3. 

where most other state have become separable. The left 
plot of Fig.[4]shows the analogous situation for decay with 
q = 1/10 and initial states with initial tangle tq = 8/3. 
Here, we have been choosing a very small value of q, since 
in this case the different excited states have comparable 
decay times so that the formation of a pointer basis is 
not very pronounced. Therefore, this example is a rather 
challenging test for our our approach that assumes the 
existence of such a basis. In contrast to the left plot, 
there is no gap between the time evolutions of the ran- 
dom states and the optimized one, and the inset that 
shows a zoom for the short time behavior indicates that 
there are a few states that have a slightly more robust 
entanglement than the state that was found ideal. How- 
ever, the differences between the tangle of the optimized 
state and the tangle of the most robust random states 
is within the error margins of the different estimates of 
mixed state tangle. 



Another feature that strikes the eye here is the fact 
that in both insets of Fig. 0] the tangle of the initially op- 
timized states turns out to be most robust over all times, 
although the optimization has only been performed for 
the initial time step. Whereas this is a mere observation 
here, a more rigorous footing the generality of such be- 
havior - that would exceed the scope of the present letter 
by far - will facilitate the search for robust mixed states 
tremendously. In particular, prior observations that the 
dynamics of entanglement 19] is often characterized very 
well by the first infinitesimal time step strongly substanti- 
ate our conjecture that optimal mixed states are actually 
given by the decay products of initially pure states. 

Here, we have been focussing on bipartite system. The 
same ideas of maximizing the temporal increment of tan- 
gle can also be applied to multipartite systems for which 
Eq. (fTJ) can be generalized [20J. The investigation for 
bipartite systems was facilitated with the Schmidt de- 
composition. For general multipartite systems, however, 
there is not a simple generalization of this tool, and there 
is an abundance of different classes of multipartite entan- 
gled states [U [22| that draws a very intransparent pic- 
ture of entangled states. Similarly to the bipartite case, 
where the Schmidt decomposition could be reproduced 
via the most robust states under very general situations 
of environment coupling, we expect that similar inves- 
tigations on multipartite systems will allow to identify 
robust classes of entangled states and provide a signif- 
icantly more physical characterization of multi-partite 
states than a mere distinction of mathematical classes 
can provide. 

Stimulating discussions with Andreas Buchleitner and 
financial support by the Alexander v. Humboldt founda- 
tion is gratefully acknowledged. 
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